Light front QEDi + i at finite temperature 
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We investigate thermodynamic properties of quantum electrodynamics in 1 + 1 dimensions 
(QEDi+i) utilizing light front dynamics. Therefore we derive the partition function of the canonical 
ensemble in discrete light cone quantization, and calculate the thermodynamical potential. This 
central quantity is evaluated for different system sizes and coupling strengths. We investigate the 
continuum limit and the thermodynamical limit and present basic thermodynamical quantities such 
as pressure, energy, and entropy, as a function of temperature for the interacting system. The results 
are compared to the ideal bosonic and fermionic cases. 
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Recently, Light Front Quantization (LFQ) originally 
introduced by Dirac in 1949 [lj has been successfully ex- 
tended to relativistic systems of finite temperatures T, 
and also chemical potentials \i @, H, 0, H, S 0, H| The 
interest in utilizing LFQ in the area of quantum statistics 
and thermo-field theory is motivated through the need to 
investigate strongly coupled relativistic systems. A perti- 
nent example of such a system is given by the quark gluon 
plasma (QGP), that is a new state of matter recently dis- 
covered at the relativistic heavy ion collider (RHIC) [9J]. 

Practical use of LFQ as a nonperturbative method 
to treat strongly interacting systems has been realized 
through Discrete Light Cone Quantization (DLCQ) in- 
troduced in [l(|. Subsequently DLCQ has been used 
to treat two-dimensional gau ge theories like Quantum 
Electrodynamics QEDi+i [ill \v\ and Quantum Chro- 
modynamics QCDi + i [lj]. Following these promising 
applications of DLCQ many other models of various di- 
mensions have been investigated. Further progress in the 
Light front Schwinger model has been obtained utilizing 
the Tamm-Dancoff approximation [3, EH ■ Besides the 
LFQ approach finite lattice size computations have been 
performed, see e.g. Refs. [H, E3 and refs. therein. Al- 
ternatively, an instant form Hamiltonian lattice approach 
in a fast moving frame was shown to give reasonable re- 
sults [HI]. The interest in QEDi + i is motivated by the 
appearance of phenomena like chiral symmetry break- 
ing, charge confinement, and topological properties asso- 
ciated to 9 vacua similar to full QCD. In addition, it is a 
widely studied theory to investigate new solution meth- 
ods like the one presented here. 

In the initial work of Ref . [2| that focused on finite tem- 
peratures, DLCQ was merely introduced to compute the 
mass spectrum of QEDi + i. Subsequently, the mass spec- 
trum has then been used in an instant form framework 
to determine thermodynamic quantities like the partition 
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function, internal energy, specific heat, etc. In this frame- 
work the authors obtained remarkable results including 
a second order phase transition at around T ~ 1 g/ \pK 
(where the strength g will be specified below). Extrapo- 
lation to the continuum limit leads to a lower bound for 
the critical exponent a > 0.7 and critical temperatures 
T c = 0.5 . . . 1.0 gj\pK (depending on the interaction) even 
in the strong coupling regime, where the Schwinger model 
is exactly solvable by bosonization. This is astonishing, 
since in a correct treatment of the two limits, (i) weak 
coupling m/g — > oo and (ii) strong coupling m/g — > 
(where m is the mass unit, see below), a free fermion 
respectively a free boson equation of state (EOS) should 
be realized. 

Within the same framework of Q Hiller et al. inves- 
tigated a super-symmetric model in two dimensions (l9l . 
|20( . They computed thermodynamic observables of the 
bosonic and fermionic sectors of the theory. It turns out 
that irrespective of th e sp ecific findings of Ref. [2| the 
models under concern [ijj |2(| due to the large N c limit 
lead to a free theory of mesonic states. As a consequence 
the EOS is trivially that of an ideal system of bosons. 

In order to study the non-ideal case we investigate the 
massive chiral Schwinger model for interaction terms of 
finite values of m/g. We first solve the DLCQ Hamil- 
tonian up to some harmonic resolution and then evalu- 
ate the proper partition function for finite temperature 
without introducing any quasi particle concept. The con- 
tinuum and the thermodynamic limits are carefully per- 
formed. 

We introduce the massive chiral Schwinger model in 
DLCQ following Refs. [H El- Within the DLCQ ap- 
proach the light-like physical extension x~ is restricted 
to —L/2 < x~ < L/2. Therefore by demanding (anti)- 
periodic boundary conditions of the fields the conju- 
gate momentum variable p + becomes discrete = 
^-n, neff. For periodic boundary conditions the n = 
mode (fermionic zero mode) exists while for anti-periodic 
boundary condition n > holds. In the following we 
choose light cone gauge A + = and neglect the fermionic 
and gauge zero modes to not obscure the thermodynamic 
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FIG. 1: Invariant mass spectrum of QEDi+i for given har- 
monic resolutions K and m/g = 1. mo denotes the lowest 
eigenvalue for K — » oo and is used for normalization. 



issues addressed here. 

The light cone momentum operators can then be writ- 
ten as 



P- 



(1) 



with the harmonic resolution K being dimensionless and 
the LF Hamiltonian H of dimension mass squared. Ex- 
pressed in Fock space creation/annihilation operators for 
particles (6^) and antiparticles [d]/) the harmonic resolu- 
tion and the free (dimensionless) Hamiltonian H are 



K=Y J n i b l b n + dldn) 
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with a high momentum cutoff A. The full Hamiltonian 
is given by [ll[ 



H = m Hq + —V = g\- 

TT 



1 



H + -V), (4) 



where m is the bare mass and g the interaction strength. 
We have introduced the ratio m/g that for m/g — > gives 
the strong and m/g — > oo the weak coupling limit. The 
Fock space representation of the potential V is a rather 
lengthy expression that has been given in Ref. ll| Sec. 
2, and will not be repeated here for the sake of brevity. 

We first determine the invariant mass spectrum M 2 = 
P + P~ = KH(K). To do so, we construct the Fock 
space for a given resolution K and collect all possible 



integer partitions {K}, where integers represent a single- 
particle momentum state of an electron or positron for 
a given momentum while the Pauli principle is observed. 
The Hamiltonian is block-diagonal in K and the result- 
ing Fock space spanned by {A"} is finite. We then cal- 
culate the DLCQ QEDi + i Hamiltonian matrix H(K), 
which can be diagonalized for values of K < 35 without 
much effort on today standard PCs. The mass spectrum 
M = \J KH(K) is shown in Figure [1] as a function of K. 

At larger K > 35 the dimension of the mass matrix 
soon becomes rather large. We therefore restrict our- 
selves to the lower part of the mass spectrum. To this 
end we introduce standard Krylov space methods to de- 
termine the lowest eigenvalues as used e.g. in Fig. [2] 

Since M 2 = KH(K) does not explicitly depend on L, 
the continuum limit is achieved by K — ► oo which in 
turn implies for any given but (due to covariance) arbi- 
trary P + > that L —t oo (continuum). Assuming some 
low order fits explained below the continuum mass spec- 
trum including an educated estimate of the error related 
to variation of the fit function or/ and parameters is read 
off at 1/K — > 0. This is demonstrated in Figure for the 
binding energy E = M— 2m of the six lowest mass eigen- 
values M for the case m/g = 1. In Table U we compare 
our results with earlier calculations [ll|, Q, [l(| Qj], HU 
that have been given for the lowest two mass eigenvalues, 
the ground (vector) state and the first excited (scalar) 



TABLE I: The different estimates for the bound state ener- 
gies Ei/g, Ev/g of the lowest two mass eigenstates at various 
m/g. For comparison the results of [lj and [TIJ using LCQ, 
finite- lattice results of [TtJ (first state), [iiil ] (second state), 
and data obtained by the fast-moving frame lattice Hamilto- 
nian method [18J are shown. Error given due to variation of 
extrapolating functions (The case g/m = 2° is also shown in 
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FIG. 2: The six lowest mass eigen values of QEDi+i as a 
function of 1/K for m/g = 1. The dashed line is a quadratic 
fit to the data (values used in the fit are colored in red) and 
to extract the continuum limit. 



state, for different interaction strengths m/g. 

We use different extrapolation functions to extract the 
low lying masses depending on the coupling regime m/g 
to arrive at our results of Table |U At low coupling (that 
is m/g > 2 5 ) the masses vary strongly and follow a linear 
curve in 1/ K, which is used in the fit, only for rather large 
resolutions K > 40. In the intermediate coupling regime 
2 _1 > m/g > 2 4 the light cone approach produces the 
most accurate results even for moderate K . However, a 
linear continuum extrapolation systematically gives too 
high estimates. Thus omitting the linear fit we do the 
analysis using second order polynomials. At small m/g 
ratios the scaling behavior is considerably different from 
the previous ones. Even at large resolution K the com- 
puted masses are still quite far away from the very pre- 
cise continuum values found in [l|| [TtJ • Changing the fit 
function to a polynomial in 1 / \f~K we obtain reasonable 
results at m/g = 2~ 2 but unstable ones at m/g = 2~ 3 . 
Therefore we do not indicate errors for the smallest m/g 
values but simply present our best fit estimates. All other 
errors have been determined through successive fits with 
next and next-to-next order power function applied to 
the same sample of data points. It is a well-known fact 
that ordinary LFQ does not capture all physics in the 
massless limit, e.g. ignores left moving fermions, that 
indicates why this approach has difficulties reaching the 
continuum. This can be cured by introducing a second 
light front [2l[ or near light cone coordinates [22]. Fur- 
thermore a proper inclusion of the gauge field zero mode 
following [231 ] should improve the convergence at small 
fermion mass. Interestingly, the bosonization of the mas- 
sive chiral LC Schwinger model as suggested in Ref. 
leads to different induced inertias, but otherwise leaves 



the interaction V in ([4]) unchanged. The bosonized model 
gives exact results in the massless limit, however, even for 
strong coupling, similar uncertainties are present in the 
fermionic approach. 

After solving the isolated DLCQ case to the present 
state of the art precision, we now turn to a canonical 
ensemble subject to the DLCQ Hamiltonian of QEDi + i 
given above. To calculate thermodynamic properties we 
evaluate the respective partition function that is the cen- 
tral quantity. Following the standard quantum statistical 
approach the partition function is evaluated by perform- 
ing a trace over the statistical operator in the discretized 
Fock space. In LFQ the respective trace is given by 

Z = Tr cxp{-^(P+ + P~)} 

= X>p{-ff X } Ck, where (g) 

{K} »• J 

The last two lines refer to the DLCQ expression, wherein 
the first part is the sum over the complete Fock space, 
and (k contains all Fock state partitions belonging to 
the given resolution K. The inverse temperature (3 = 
l/(ksT) (Jcb Boltzmann constant) is the Lorentz invari- 
ant rest-frame temperature. In order to keep the stan- 
dard definition of the temperature the (however kine- 
matic) P + component in the proper definition of the sta- 
tistical operator on the light front appears in addition of 
the light front Hamiltonian P~ . Note that only for the 
interaction free case (P + + P~)/2 = P° can be replaced 

by P° = V P 2 + M 2 , where the Lorentz invariant mass 
eigenvalues M may be evaluated is any of the relativis- 
tic forms, e.g. light front form, as has been utilized in 
Refs. @, [ljj[2(|. The case considered here, however, in- 
cludes interactions and hence the evaluation of the trace 
becomes more cumbersome. 

To give an instructive example, which is also used to 
estimate the theoretical error of the calculation, we evalu- 
ate the thermodynamic potential density u>f = Clf /L for 
the canonical ensemble of the ideal fermi gas. The ther- 
modynamical potential for two particle species is given 
as 

oo 


(6) 

where we set the mass to m = 0.5MeV. To compare © 
to the standard instant form expression, the r.h.s. can 
be directly transformed by a simple variable substitution. 
In the thermodynamic limit, the instant and light front 
form give equal results for the ideal case. 

Evaluating the free thermodynamic potential density 
(J5]) for a given temperature will be referred to as the an- 
alytical result. To demonstrate our method this will be 
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FIG. 3: The thermodynamical potential f2/ = — 2irTlnZ as a 
function of L. Open symbols depict the potential at maximal 
resolution K = 110 for the different temperatures. Closed 
symbols are given by the extrapolation as explained in the 
text. The slope of the linear part (values selected colored in 
red) is fitted to extract the invariant potential density. 



compared in the following to the direct numerical evalu- 
ation of (JSJ) for V = that depends on the length scale L 
and a maximum K. Since In Z is an extensive property 
it scales with L as well as ui = —T\nZ. For a given max- 
imal harmonic resolution K the numerical results leads 
to InZ as a function of the system size L. Hence InZ 
can be read off the slope with respect to L. The situation 
for different temperatures is depicted in Figure [3] The 
analytic (exact) results are given as dashed lines and the 
numerical values as symbols, empty ones reflecting raw 
numerical data up to K = 110. The filled symbols are 
raw data corrected by a nonlinear fit as explained below. 
For large system sizes the harmonic resolution used in 
the evaluation is not sufficient, hence a deviation from 
the expected linear dependence of InZ on L appears. 
For small L the system is too small to read off the ther- 
modynamic limit, which also leads to a derivation (finite 
size effects). There is an optimal region of system sizes 



TABLE II: Values of In Z for the free fermi gas extracted from 
a linear fit to the numerical analysis along with the analytical 
results and the estimated uncertainties. 
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-0.862941 


-0.863176 


0.03 


0.9 


-2.24692 


-2.24824 


0.05 


1.2 


-4.24937 


-4.24566 


0.07 


1.5 


-6.85920 


-6.76379 


1.39 



where InZ depends almost linear on L, which is used to 
extract the slope by making a linear fit (indicated by red 
symbols in Fig. [3]). In order to improve the numerical 
results we have used the following fitting algorithm for Z 
of ([5]) : Beyond the maximum of Z we approximate C,k by 
an exponential of ^2 n= _2 a nK n . The function depends 
on five parameters a n that are determined by a x 2_n t to 
the largest 30 calculated K-values and then used for an 
extrapolation to larger K. The resulting values of InZ 
for the four temperatures are given in Table HT1 The error 
indicated is the difference of the numerical from the ana- 
lytic result. For small L the analytic value is approached 
from above and for large values it deviates again to larger 
values. As a consequence the slope (i.e. InZ) is usually 
slightly too steep. For the interacting case, where no 
exact analytic result is available, we utilize the same fit- 
ting procedure. Note however two differences to the pure 
fermionic case. First, particle number is not conserved 
because of an interacting field theory and, second, the 
charge Q of each state of the statistical ensemble is set 
to zero due to confinement. To be more explicit, the ideal 
fermionic (bosonic) gas as a limit of g/m — > (m/g — ► 0) 
has been evaluated for chemical potential [i — 0. In the 
fermionic case this has to be done using the full numer- 
ical calculation obeying the condition Q = 0. However, 
the bosonic case can be worked out directly by evaluating 
the following integral 

oo 


(7) 

since Bose-Einstein condensation is absent in one dimen- 
sion [24| (and would occur only at fi = m for an ideal, 
relativistic bosonic gas in three dimensions). 

In the reminder we give several thermodynamic quan- 
tities. Note that volume and temperature are now given 
in units of the ground state mass mo that sets the scale 
of the interacting system. In Fig. 2] we show the pres- 
sure p = —lu = TlnZ/L plotted as a function of tem- 
perature. Division is by T 2 (rather than the usual T ) 
is due to the reduced dimension of the problem. The 
line shows the analytic result for the ideal bosonic case 
(limit m/g — ► 0) according to |(7J), the diamonds the ideal 
fermionic case (limit m/g — > oo). The other symbols rep- 
resent two different interacting cases. The strong inter- 
acting one m/g = 2~ 3 (squares) is closer to the bosonic 
case, whereas the weak interacting one m/g — 1 (circles) 
is closer to the ideal fermionic case. There is a smooth 
transition from the bosonic case to the fermionic case, 
depending on the interaction strength. In particular, we 
do not observe a critical behavior seen in [3] that might 
have been due to the rather low harmonic resolution used 
there. For T — > oo, i.e. the classical ultrarelativistic case, 
masses can be neglected m — > and hence all curves ap- 
proach p/T 2 = 7r/6. From Table U we assume a conserva- 
tive mean error of 2.5% for the interacting case indicated 
in figures by error bars. 
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FIG. 4: Pressure as a function of temperature divided by T 
for different models: pure bose gas (solid line), stronger inter- 
acting Fermi system m/g = 2 -3 (squares), weaker interacting 
Fermi system m/g = 1 (circles), noninteracting Fermi gas in 
Q — sector as explained in the text (diamonds). Estimated 
error as explained in the text. 




FIG. 5: Internal energy density of interacting gas as a function 
of temperature divided by T 2 , for explanation of curves see 
caption of Fig. [4] 



The energy density u is achieved by thermodynamic 
relations involving a numerical derivative of In Z 1 hence 
error bars indicated are bigger. With the same coding 
u/T 2 is shown in Fig. In the classic ultrarelativistic 
limit the energy density converges to u/T 2 — > w/6. 

The entropy density s = (p + u)/T divided by T is 
shown in Fig. [5] with the classical ultrarelativistic limit 

s/r-> tt/3. 



FIG. 6: Entropy density of interacting g cLS ELS £L function of 
temperature divided by T, for explanation of curves see cap- 
tion of Fig. m 
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FIG. 7: Equation of state, pressure p vs. energy density u for 
QEDi+i, for explanation of curves see caption of Fig. [4] The 
data points belong to the following temperatures T/mo = 
0.3,0.6,0.9,1.2,1.5. 



The resulting equation of state is shown in Fig. [7] For 
moderate temperatures p ~ 2u, which is taken to esti- 
mate the error in u to be twice that of p. 

In conclusion, we have presented for the first time a 
possible way to directly evaluate the partition function 
of a strongly interacting relativistic system within the 
framework of light front quantization. To this end we 
have utilized DLCQ that leads to a Hamiltonian evalu- 
ated in a Fock space of finite large dimension. As an 
interaction we have chosen QEDi+i. For illustration 
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we have investigated some intermediate cases between 
the noninteracting (fermionic) limit (chemical potential 
fx = 0, charge Q = sector) and the completely coupled 
(bosonic) limit of QEDi + i. 

Present numerical effort is manageable, i.e. all calcula- 
tions have been performed on a regular PC. It is also de- 
sirable that besides evidence of feasibility further thermo- 
dynamic quantities could be calculated, such as specific 
heat or speed of sound. Some additional effort is needed 
here, since either for use of thermodynamic relations (nu- 
merical) derivatives are involved, or some quantities like 
energy U — (H) have to be calculated from thermody- 
namic averaging. In any case within this framework it is 
also possible to evaluate correlation functions and deter- 



mine thermal masses in a next step. 

In view of demands to investigate the QCD phase tran- 
sitions, the method can now be extended to QCDi + i, 
e.g., along the lines of [13j. A further next step along 
this line would be the inclusion of transverse degrees of 
freedom, and finally the extension to finite chemical po- 
tentials which is by construction not biased by fermion 
doubling and sign problems due to the generic Hamilto- 
nian approach. 
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